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Abstract 

A consistent census of penguins in the AI = 1/2 rule is taken from the tjq pole contribution 
to the radiative Kl — ► 77, Kg ~~ > vr°77 and K + — > 7r + 77 decay modes. We briefly comment on 
its impact for Kl — ► 7r°7r°77, Kl — ► tt + tt~j and check its compatibility with the Kl — Kg mass 
difference and the CP violating e'/e parameter. 
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1 Introduction 



A precise quantitative understanding of hadronic kaon decays has, up to now, been upset by the 
non-perturbative nature of strong interactions at low energy Though qualitatively, there is little 
doubt that the AI =1/2 rule is a pure QCD effect, the genuine mechanism still eludes us. In 
particular, the relative strength of penguin and current-current operators has been debated for 
years. This has in turn impeded theorists from making use of the remarkable precision achieved 
in e'/e measurements. 

In the present work, this question is addressed from a phenomenological point of view. The 
relative contribution of penguin and current-current operators in the AI =1/2 rule is known 
to be accessible from anomaly-driven radiative K decays, i.e. decays occurring through ir, r/, r( 
poles. Still, a thorough investigation in the context of large-./V c Chiral Perturbation Theory 
(ChPT) is called for. In particular, we will see that U(3) ChPT holds the key to understand, 
and go beyond, the well-known vanishing of the SU(3) Kl — > 77 amplitude at lowest order. 
Besides, the experimental situation has been improved and is expected to continue doing so for 
a number of radiative K decays. It is thus appropriate to review their theoretical treatment, 
especially regarding r/o effects, and to guide experimentalists towards suitable observables to 
probe the AI = 1/2 rule. 

The paper is organized as follows. In the first section, the large- N c formalism is recalled, 
with emphasis on the U(3) ChPT weak operator basis. In particular, the strategy to extract the 
penguin fraction from phenomenology is exposed. The next few sections deal with radiative K 
decays. We first review in great details the Kl — > 77 mode, and develop the phenomenological 
tools to deal with 7/ — rj mixing to be used in the rest of the paper. Then we go on with the 
Ks — ► vr°77 decay and the pole contribution to the K + — > 7r + 77 one, which are both quite 
sensitive to 77 — rf effects. The Kl — > 7r°7r°77 and Kl — ► vr + 7r _ 7 modes are only briefly analyzed 
as no new information can be obtained from them. Finally, the implications for the hadronic 
observables AMls and e'/e are investigated. Our results are summarized in the conclusion. 

In the appendix, the complete analysis of Kl — ► 77 in SU(3) at 0{p G ) is presented, some 
aspects of which are used throughout the paper to investigate the reduction of U(3) amplitudes 
to SU(3) ones. 



2 Theoretical framework 



Our starting point is the QCD-induced AS = 1 effective weak Hamiltonian below the charm 
mass scale [1,2]: 



nfff 1 (/x < m c ) = ^v ud v: s £ 

i=i 



Qi 



v ud v* s 

with the familiar current-current (i = 1,2) and penguin (i = 3,..., 6) operators 



Qi = 4 (sLj a d L ) (u L j a u L ) , Q2 
Q3 = 4 (sLj a d L ) {qLl a qL) , Qi 
Q5 = 4 {sL7 a d L ) (qRl a qR) , Qa 



^{sl1o.ul) {uLl a d L ) , 
4(s L 7 Q g L ) {qLl a dh) , 
-8(s L qR) {qRd L ) , 



after Fierz reorderings. In our notations, q^ = |(1 ± 75)^ and the light flavors q = 
are summed over. The connection between Eq.Q and kaon phenomenology requires of 



u, d, s 
course 
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additional tools, QCD perturbation theory being helpless to estimate the Qi hadronic matrix 
elements. We will make use of chiral Lagrangian techniques, as we now describe. 



2.1 U(3) chiral representation of AS = 1 weak operators 

As is well-known, the above operators are not all independent. Indeed, we have the identity 

Q 2 - Qi = Qa - Qs- (3) 

Besides, Q4 and Qq have the same color and flavor structures. The chiral realization of the 
effective Hamiltonian 7~L^ff =1 must thus allow for an explicit representation of the left-handed 
flavor-singlet current in Q3 if we want to be in a position to disentangle Qip and Qe contributions 
to the AI =1/2 kaon decay amplitudes. For this reason, we will perform the hadronization of 
Eq.Q in the U(3) chiral expansion, with an explicit flavor-singlet degree of freedom r/o, rather 
than in SU (3). 

The extension of SU(3) to U(3) proceeds through the large-iV c limit, with N c the number 
of QCD colors. Considering iV c as large is the key to a consistent disposal of the QCD U(1)a 
anomaly [3]. The spontaneous symmetry breaking U(3)l x U(3)r — > U(3)v gives then rise to a 
nonet of pseudoscalar Goldstone bosons, which are written 
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in the exponential parametrization. The anomalous breaking of the U(V)a symmetry is reintro- 
duced through 1/N C corrections. In particular, the leading nonlinear Lagrangian reads [4] 



C 



(p 2 ,oo)+(p°,l/Af c ) 



^{d^ud^) + ^{ x rf + u x ] ) 



+ 



16iV, 



■mg(ln[/-ln C/ f ) 2 



(5) 



where () denotes a trace over flavors, the external source x is frozen at x = r M to account for 
meson masses, M = diag(m u ,md,Tn s ) is the light quark mass matrix and mo represents the 
anomalous part of the r/o mass. The ChPT symmetry-breaking scale F is identified with the 
neutral pion decay constant F n = 92.4 MeV at this order. Note that the leading SU(3) chiral 
Lagrangian is recovered in the limit mo — ► 00, when the r/o decouples. 

The effective Hamiltonian (Q) contains both (8z,,1r) and (27l,1r) representations of the 
chiral group U(3)l x U{3)r. At leading order in the momentum expansion, four operators built 
out of U and X have either of these transformation properties and at most one factorized r/o 
field 1 : 

Q 8 = 4 (L^) 23 , Q| = 4 (L„) 23 (L») , Qf = F 4 ( x E/t + urf) , 

(6) 

Q27 = 4 [(L„) a3 + § (L,) 13 m 2l - \ (L M ) 23 (L»)] , 

where the matrix collects the chiral realizations of the left-handed currents from Eq.©: 



i—ld^UuA =(Lj k (q l > k =u,d,s). 



(7) 



1 Operators with two or more rjo factors have negligible effects on the processes we will consider due to isospin, 
loop and/or 1/N C suppressions. 
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Qs, Q™ and Q27 are the trivial U(3) generalizations of the usual SU(3) octet and 27 operators 
while 

Ql ~ (L M ) 23 3% (8) 

is peculiar to the U(3) framework. Note that the introduction of one extra (InU — InU^) ~ 770 
factor in Qs, Q™ or Q27 would violate CPS invariance. Besides, other structures like (xU^ — 
^X t )2 3 (ln?7-lnC/t) can be brought back to Q$ through the use of the classical equations of 
motion. For our purpose, the complete low-energy realization of TL^p 1 at 0(p 2 ) is thus given 
by 

7i w ~ = G 8 Qs + G S 8 Q S 8 + G™Q™ + GiiQtj- (9) 

Yet the weak mass term Q™ is known not to contribute to hadronic observables at leading order, 
while its higher order effects can be absorbed into a redefinition of the weak counterterms [5,6]. 
We will thus disregard it from now. Later, we will check explicitly that the same holds true in 
respect of radiative kaon decays. 

Let us now make a closer connection with the current-current and penguin operators of 
Eq.@. 



2.2 QCD-inspired alternative basis 

As a first step towards an alternative QCD-inspired formulation of Eq.Q, we consider the 
effective Hamiltonian (^Q) in the factorization approximation. 

In that limit, the chiral realization Qi of the (V— A)<g>(V±A) Qi operators is straightforward. 
It simply follows from the hadronization of the left- and right-handed currents in Eq.0 (with 
U <-> U' for the right-handed currents). This gives: 

Qi = 4 (L M ) 23 (L») u , Q 2 = 4 (L M ) 13 (L») 21 , 

Qs = ~Q 5 = 4: (L„) 23 (LP) , g 4 = 4 (l^) 23 . 

For the Q§ penguin operator, the situation is more delicate as the constraint UW = 1 forces us 
to go beyond Eq.© to hadronize the quark densities [7,8]. From the 0(p 4 ) strong Lagrangian 
in the large- N c limit 

4 p4,oo) 3 (d»ud»u\ x u ] + u x f )), ( x arf + au x ^), (( x u ] ± u x f ) 2 ), (n) 

two structures 

Q 6 = 4(L^) 23 (12) 

(identical with Q4) and Qg 1 emerge. 

Corrections with respect to the factorization approximation consist in meson exchanges be- 
tween currents or densities. Such 1/N C loop effects, despite formally suppressed, are important 
due to the quadratic dependence on the physical cut-off for a truncated chiral Lagrangian de- 
scribing the low-energy strong interactions of massless pseudoscalars in terms of the scale F [9] . 
Accordingly, no new structure will be induced at 0(p 2 ) since the other scale parameter ijiq in 
Eq.© has no effect on this fast operator evolution [10]. We thus reach the ansatz: 

6 

Hfa s=1 = GwY,*iQi- (13) 

i=l 
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The overall constant Gw and the weights Xi are not fixed by symmetry arguments. Let 
us set Gw = GfV u< iV* s /\/2 = 1.77 x 10~ 12 MeV~ 2 so that in the factorization approximation 
(neglecting V td V t *): 

FACT. , > C,a\ FACT. 2 F K / F^ - 1 

Xi = Zi{Hf act ) ijfcO, x 6 = -r — 2 T z 6 Wact) > ( 14 ) 

m K — m% 

with Hf ac t around 1 GeV. The difference Fx/F n — 1 is generated by a linear combination of the 
first two terms of Eq. (|lH) . 

The actual Xi values differ from those given in Ea. (|14|) by long-distance strong interaction 
effects. Following the principles of chiral expansions, we wish to extract them from phenomenol- 
ogy. Yet only three combinations of x% are accessible as the corresponding operators can all be 
expressed in terms of Qg, Q§ and Q 2 7- Identifying Eqs.Q and (|T3|) . we obtain: 

2 3 

Gg/Gw = - -=xi + -x 2 +x 4 + x 6 
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G s 8 /G w = -^xi - -x 2 + x 3 - x 5 (15) 
3 

GtjIGw = -=(xi + x 2 ) ■ 
5 

In order to invert the above system, some extra theoretical assumptions are needed. Corrections 
with respect to the factorization limit are likely to be important in the case of xq due to the 
presence of the factor r = 2m 2 K /(m s + m^) in Eq. (j!4() . We thus have to keep it. However, the 
situation is different for x$, x^ and x$, whose small factorized values suggest marginal effects. 
Neglecting Q3, Q4 and Q5, the QCD-inspired set (Qi,Q2,Qe) becomes equivalent to the basis 
(Qs; Qgi Qm)- Notice again that, had we worked in the SU(3) framework, the disappearance of 
the second octet would have led to a linear relation between Qi, Q2 and Qq, rendering thus 
impossible the distinction between current-current and penguin operators from phenomenology. 

2.3 A first look at the weak couplings 

Our goal is to extract x\, x 2 and xq. A first piece of information can readily be obtained 
from the analysis of the K — > irir decays. In the isospin limit (m u = rrid, no electromagnetic 
corrections), these can be parametrized in terms of two isospin amplitudes, Aq (AI = 1/2) and 
A 2 (AI = 3/2): 

A (kO ^ ^-) = Ao e iS + -~^A 2 , 

\/2 

A (K° -> 7r°vr ) = A e iS - V2A 2 , (16) 

a(k+^^^) = \a 2 , 

with 5 = 5o — 82 — vr/4, the final state interaction phase shift. In our conventions, A2 is real and 
positive in the limit of CP conservation, adopted from now on. The AI =1/2 piece receives 
contributions from both the Qg and Q27 operators (none from Q| in the isospin limit), while 
the AI = 3/2 piece is only brought about by Q27- In terms of the Xj's, this gives: 

A) = V2F (g 8 + ^G 2 ^j (m\ - ml) = ^-FG W (- Xl + 2x 2 + 3x 6 ) (m 2 K - ml), (17) 

10 2 

A 2 = —F G 27 (ra\ - ml) = -F G w (x t + x 2 ) (m 2 K - m 2 ), (18) 
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the QCD-inspired estimate of Eq. ([23|) , 



with X4 neglected (as indicated above) and no contribution from 2^5. Experimentally, uj^ 1 = 
^0/^-2 = 22.2, which is the very statement of the AI = 1/2 rule. This large factor implies a 
dominance of Gg over G27: 

G 8 = 9.1 x 1(T 12 MeV~ 2 , G 27 = 5.3 x 1(T 13 MeV~ 2 . (19) 

However, this does not give any information yet on the relative sizes of the x^s, and thus on the 
current-current and penguin fractions in Aq 

Tec = ~ X1 + 2X2 and = ^ . (20) 

—x\ + 2x 2 + 3x 6 —x\ + 2x 2 + 3x 6 

The missing piece of information is of course G|. Before proceeding to its phenomenological 
determination, let us see what comes out from additional QCD-inspired assumptions. The values 
of xi, X2 and xq corresponding to a given ratio G^/Gg after imposition of the constraints (|T9|) 
on the system (|T3|l are depicted in Fig.l. In the case of x\ and X2, strong interaction effects are 
not expected to change the signs of the factorized values qi4|) . i.e., we should keep x\ < and 
X2 > 0. This favours a negative ratio G|/Gs, excluding large positive values of xq. A natural 
scaling further points towards | C| / Cg I < 1> which would guarantee Xi ~ 0(1) and also preserve 
the sign of x% with respect to the factorization approximation, i.e., xq > 0. A small fraction 
of penguins, as in the factorization approximation at 1 GeV where Tp is around 20%, would 
require xq ~ 1, i.e., a ratio G^/Gg as negative as —0.8. However, another theoretical clue can 
be obtained from the nonlinear relation between the Wilson coefficients zi(fi) and Z2 (fi) derived 
in the leading logarithmic approximation [1]: 

(z 2 + Zl ) 2 (z 2 - Zl ) = 1. (21) 

Indeed, this relation is known to receive rather small corrections at next-to- leading order (even 
for [x as low as 700 MeV). Assuming 

(x 2 + Xl f (x 2 - Xl ) = 25Gh f 7 2 J 3 ~ 6GI) - 1 (22) 
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to within, say, 30%, would then give the smaller (in absolute value) ratio: 



(G s 8 /G 8 ) th = -0.38 ±0.12. (23) 
2.4 Tracking down penguins at the poles 

Let us now tackle the phenomenological determination of G|. The vertices induced by Q| always 
involve one single rjo field (see Eq.(jHJ). It is thus natural to turn to r/o pole contributions to 
anomaly-driven radiative kaon decays. These proceed through the 0(p 4 ) Wess-Zumino-Witten 
(WZW) action of U(3) ChPT, which is for QED external sources [11]: 

£wzw = ^2 e e^A^d v Ud p ^d a \J{U\ Q}) , (24) 

4W = ^2^ 2 e^F^Ap ((QQ{d a U, l7+» + \{QU^Qd a U - QUQd a U^ , (25) 

with A p , the electromagnetic field, = d p A u — d u A p the corresponding strength tensor, 
Q = diag(Q u , Qd, Q s ), the light quark charge matrix in units of the positron charge e and s^ pa , 
the Levi-Civita tensor with e 0123 = +1. Note that in f7(3) ChPT, there are additional unnatural 
parity operators at 0(p 4 ) (necessarily chiral invariant as unrelated to the WZW action), like for 
example 

e^F^Ap (QQ) (d a UW) , (26) 

but, being 1/N C suppressed, they will be discarded. This is supported by the analysis of the 
decays r/, rf — ► 77, reasonably well reproduced at leading order by Eqs.(|25j) and ® alone [12]. 

In the next sections, we will consider pole-dominated radiative modes like Kl — > 77, K$ — > 
7r°77 and Kl — * 7r°7r°77, as well as transitions that receive also other types of contributions like 
K + — > 7r + 77 and Kl — > ir + ir~j. 



3 The Kl —> 77 decay 

Despite a long history, and precise measurements by the NA48 and KLOE Collaborations, 
T {K L -> 77) / (2.81 ± 0M stat ± 0.02 syst ) x lO" 3 [13] 



T (K L -> vrOvrOvrO) \ (2.79 ± 0.02 stat ± 0.02 syst ) x lO" 3 [14] ' (27) 

it is fair to say that a satisfactory theoretical description of the decay Kl — * 77 is still lacking. 

At lowest order in the chiral expansion, i.e. 0(p 4 ), it proceeds through pseudoscalar poles, 
as depicted in Fig. 2a. In the SU(3) framework, the rjo is a higher order effect and only the tt° 
and 778 are allowed to propagate, with the result that the amplitude vanishes exactly: 

A» ( Kl - 77) SU = ] — {G27 ~ G s + G™) m\ ( 2 1 + 2 1/3 2 ) is^k lp k 2a = , 

7r V m K — m~ m K — m* 8 I 

(28) 

upon enforcing the Gell-Mann-Okubo (GMO) mass relation (valid at this order): 



4m^ — ml 



= -_H1K_^ _ (29) 



The SU(3) decay amplitude therefore starts at 0(p 6 ), with in particular all the effects of the 
r/o- At this stage, phenomenological pole models are constructed [15-19], involving the physical 
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rf'' masses but also chiral symmetry breaking parameters for the weak transitions Kl — ► 7r°, rj, rf 
and the radiative decays 7r°, 77, 77' — > 77. Yet no definite prediction can be attained because, 
having a vanishing lowest order, the rate exhibits a high sensitivity to these phenomenological 
parameters. 

3.1 Pole amplitude in U(3) ChPT 

In the present work, large- N c ChPT is used to introduce the effects of the rj state consistently 
in both the weak and strong Lagrangians Eqs.©, Q and l|25|). The amplitude is then non- 
vanishing already at 0(p 4: ), featuring the parameter of interest G|: 

ST {KL 77) m ^ ( G| + | G27 ) ( m§ _3^ +3m? ) . ,30, 

Yet the contributions of the Qs and Qg 1 operators cancel again once the 77s — 770 theoretical mass 
matrix is inserted 



/ 4774-77-4 2x/2, 2 , \ 



TTi 2 

mm 



3 



m; - 77%; 



,2 1 m 2 



\ -^-( m i- m K) m o + — % — / 



(31) 



which corresponds to the generalization of GMO to U(3). 

Looking back at Ea. ((T3|) . G| + §^27 is precisely the weight of the Q\ operator. This is 
the only one able to generate a uu pair from the incoming Kl at lowest order {Qz^ are again 
expected to give small contributions, and are discarded). The flavor and color structures of Q2 
prevent Q2 from contributing, while Qq appears only through its dd + ss component and cancels 
out 2 . The decay Kl — > 77 thus proceeds entirely through the long-distance sd(ds) — » uu — > 77 
transition at lowest order (Fig. 2b), which is in agreement with a quark-level analysis [16]. 

Let us stress once more why working in ?7(3) is essential. Beside the fact that the quark 
basis uu,dd,ss is ambiguous in SU(3), there are not enough \AS\ = 1 operators to permit a 
definite identification of the underlying transitions in this symmetry limit. In fact, Qi,Q2 and 
Qq are not independent in SU(3), hence the total Kl — > 77 amplitude has to vanish at 0(p 4 ) 
since the last two do not contribute. 

The above interpretation can be made obvious by a direct computation in the quark basis: 

(K L - 77) U = vLak ■ Plo (m 2 K ) ■ V% , (32) 



2 Such a cancellation of the Qs penguin operator has already been noticed in the mo — > limit [20] 
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with the LO propagator (R = 2N c (m 2 K — m 2 )/m 



o 



iPLO If)- 1 = ( 1 ) 1 1 ! | , (33) 





the WZW 77 vertices (proportional to Q 2 q ): 

V2N c a 




Vg = - :L ^~ I 1/9 I ie^k lp k 2a , (34) 



and the weak vertices, in both the (Qg, Q$, Q27) and (Qi, Q2, Q§) bases (still keeping Q 

( G% + 2G 27 /3 \ / Xl \ 

V weak = i2V2F 2 m 2 K G 8 - G^ + G§ - G 27 /3 = i2^2F 2 m 2 K G w x 6 . (35) 
\ G 8 - G™ + G| - G 27 /3 / V x 6 / 

The cid and ss components of V wea k are equal due to CPS invariance. One can easily check that 
only the uu component of the weak vector contributes to the amplitude, as announced. 

A priori, one could think that NLO effects may reintroduce significant Qq contributions. 
This intuition is however based on the pathological SU (3) situation, in which there is an ex- 
treme sensitivity due to the inability of SU (3) to catch the leading order sd(ds) — ► uu — ► 77 
contribution. Once within U(3), this piece is correctly identified and we can consider that NLO 
effects will behave according to the usual chiral counting, i.e. are at most of 30%. In other 
words, we will work under the assumption that Q\ represents the dominant contribution to the 
decay. 



3.2 Physical mass prescription 

Diagonalizing Eq. (|31j) . the leading order 77 and rj theoretical (or Lagrangian) masses are found 
to be 

TYl 2 1 / 4 

m lw) = ^+ m ^ ± 2V m °~ 3^ ~ m ^ m ° + 4 ( m ^ ~ m ^ 2 ' ( 36 ) 

The parameter tuq has to be fitted to reproduce both m T] and as close as possible to their 
physical values, which gives mo between 800 and 900 MeV. From Fig. 3, the masses are then 
well reproduced, within less than 15%, which is quite consistent with chiral counting for higher 
order effects. 

Still, in the present work, we are mainly concerned by pole amplitudes, i.e. amplitudes 
involving virtual 77 — 77' exchanges. For such processes, a good physical mass prescription must 
ensure correct analytical properties, i.e., freeze the poles at their right places. Varying tuq is 
therefore not satisfactory. In analogy with the usual 0(p 2 ) substitutions 

ml = T - (m u + m d ) -» Ml ~ (135 MeV) 2 , 

m 2 K = \ (m s + m u4 ) - M 2 K ~ (495 MeV) 2 , (37) 
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our prescription is to substitute for the LO mass matrix Ea.pi|) the most general mass matrix 
compatible with isospin symmetry 3 . In the uu, dd, ss basis considered in Eq. (|33|) . the (tt°, rjs, rjo) 
propagator takes then the generic form 

iP ph ys{ q 2 )-l = {e-Ml)[l I ]-^°( 1 1 1-8 ), (38) 

with Mq,5,R, some parameters to be adjusted to freeze the poles at their physical values. 
Equivalently, working in the rfei^o basis, this propagator can be written 

, 2 \-i _ / cos6» P sin Op \ ( q 2 - M 2 W cos<9 P -sin6p\ . . 
phys [Q UfK ~ V -sintfp cos^p J V q 2 -M 2 )\ sin0 P cos#p ) ' ^ 

and the 7r° does not mix (for the model-independent relationship between Eqs. (|38|) and (|39|) . see 
Ref. [21]). The masses are then fixed to 

M v ~ 547.8 MeV, Afy ~ 957.8 MeV, (40) 

with the mixing angle around Op ~ —22° if a large-iV c limit is taken at each order in the 
momentum expansion [21]. 

When the leading order is non-zero, inserting physical masses amounts to including a par- 
ticular class of higher order effects. Yet there are still many other sources of corrections, and to 
get a handle on their order of magnitude, we will consider 6p as an effective angle, allowed to 
vary around —22°. Specifically, the sensitivity of the pole amplitudes will be probed by taking 
6 P inside the range [-15°, -25°]. 

For example, we could include, in addition to the mass corrections, the deviations from 
F v = F„i = F n by introducing the renormalized decay constant matrix as 

Pphys (q% svo = F 2 (F,)" 1 . ( Pphys (q 2 )^) . (F^)" 1 . (41) 
This amounts to switching to a two-mixing angle formalism since by definition 
cos 6p —sin9p \ _ f F§ Fq \ _ / F$cos9g — i*osin#o 



sin'V rosfV J Fv ~\F^ tff ) ~ \ F 8 sin 6 8 F o cos0 o '' ' L> ) 



3 In principle, the pseudoscalar widths have to be included to really ensure correct analytical properties. How- 
ever, these are numerically not relevant and will be neglected. 
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However, not much is gained by using this prescription instead of Eq . (j39j) as there are still 
unaccounted NLO effects beyond masses and decay constants. Numerically, we have checked 
that in all the cases considered, taking Ea. (|42l) with the values of F%, Fq, 6$ and Oq from either 
the large N c analysis of Ref. [12] or the phenomenological fit of Ref. [22] does not alter our 
conclusions (in general, the deviation with respect to Eq. (|39j) with Op ~ —22° can be accounted 
for by taking Op ~ —15° or slightly smaller). 

3.3 The Kl — > 77 rate and its implications for the AI =1/2 rule 

Naively enforcing the physical mass prescription, i.e. taking Eq. (j32j) with the propag ator En.(l38l). 
would however reintroduce a large Qe contribution. This is a spurious effect because it corre- 
sponds to the partial inclusion of higher order terms (only those from mass corrections) for a 
contribution which vanishes at lowest order. As said before, the decay should be dominated by 
the uu transition, hence large cancellations are expected to occur for Qq at NLO (indications of 
these are discussed in the next subsection) . 

On the contrary, since the uu part is non-vanishing at lowest order, inserting the physical 
masses instead of the Lagrangian ones is allowed. Therefore, we definitively discard the dd and 
ss components from the weak vertex: 



V weak = i2V2F 2 m 2 K G w = i2^2F 




m 



K 





V 



(43) 



and use Ea.()38[) to reach the final form for the amplitude (eg = cosOp, sg = sin#p): 
A w (K L - 77 ) = ^ [g% + 2 -G 2 ^ M 2 K ie^k lp k 2a 



1 (eg - 2V2s e ){c - V2s e ) (sg + 2V2c e ){s e + V2c e ) \ 



\Ml-M 2 3(M 2 K -M 2 ) ?,(M 2 k-M 2 ) ) 

Numerically, this amplitude is dominated by the r/ pole: 

A" ( Kl 77) = (g s s + ^G 27 ^j [(0.46) w - (1.83 ± 0.30)^ - (0.12 ± 0.02) v ] ie^k lp k 2a , 

(45) 

where the numbers are in MeV and the errors amount to varying Op between —15° and —25°. 
The parametrization is rather stable with respect to rj — rj' mixing parameters, which means 
that the amplitude is under control once the poles are frozen at their right places. Then, from 
the experimental values Eqs. (|27ft with T(Kl — ► from [23,24], the only free parameter 

can be extracted up to a two-fold ambiguity: 

(G s 8 /G 8 ) ph * ±1/3 . (46) 

From QCD-inspired arguments, the positive solution is discarded as it would lead to x±, x 2 
of the wrong signs and xq very large (see Fig.l). The negative solution 

(G s 8 /G 8 ) ph+th = -0.30 ±0.05 (47) 

is detailed in Tableland Fig. 4. Also, (x± + x 2 ) 2 (x 2 — x{) stays remarkably close to one, which 
means that the low energy evolution of current x current operators is indeed rather smooth. 
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Table 1: The coupling G| as extracted from Kl — ► 77, in terms of Op, and the corresponding 
weights Xi for the current x current and penguin operators. The last two columns give the 
current x current and penguin fractions in Aq defined in Eg. 1)20(1 . 
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Figure 4: Graphical representation of Table 1: a) Weak operator weights X{ as a function of Op, 
b) Corresponding current x current and penguin fractions in Aq. 



On the other hand, xq increases quite dramatically from its perturbative QCD value Ea, (|14|) 
and appears responsible for about 2/3 of the AI = 1/2 enhancement 

T P = (68 ± 6) % (48) 

where the error reflects only the variation of Op. Let us stress again that x%, X2 and x§ are 
the weights of the weak operators at a very low hadronic scale, around the ChPT scale F. 
It is intuitively clear then that xq, and thus Tp, can be large since we allow the evolution 
Qi, Q2 — >■ Qe all the way down to this scale. 

3.4 Effective suppression of Q e at higher order 

In the context of SU(3), there is only one weak octet operator, hence it is clear that large C(p 6 ) 
"Qs" effects are present. Indeed, it is through the counterterms that the leading sd(ds) — > 
uu —* 77 transition is reconstructed. Yet, this reconstruction is not trivial and implies a specific 
interplay between strong and weak counterterms of both 0(p 4 ) and 0(p e ), corresponding to 
corrections to the decay constants, two-photon vertices, masses (i.e., corrections to GMO) and 
weak transitions. 

In the appendix, the complete analysis of Kl — > 77 at 0(p 6 ) in SU(3) ChPT is presented. 
It is shown that, to a large extent, the counterterms either act together to reconstruct Qi, or 
cancel among themselves. For example, corrections brought by F^g 7^ F n cancel out with some 
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mass corrections, and are thus at least of 0(p 8 ). Only a small irreducible combination of (scalar- 
saturated) counterterms survives, on which not much can be said at present. Still, given the 
dynamical insight gained, it is reasonable to expect this combination to be small. In conclusion, 
there is a clear indication that no large contribution from Qq is generated at 0(p e ) in SU(3) 
ChPT. 

Something very similar should happen in U(3), but it is unfortunately difficult to be as 
quantitative since we do not have any handle on the U(3) weak counterterms. Anyway, here 
also, it is likely that corrections to the masses, decay constants, weak and WZW vertices cancel 
among themselves, at least in part, for the Qe contribution. 

This is to be contrasted to the phenomenological pole model analyses [17,18], which overlook 
cancellations among 0(p & ) corrections. To be more specific, these models usually introduce a 
parameter p to account for nonet symmetry breaking. In U (3), this breaking is entirely due to the 
Qq and Q27 operators (which induce only Kl — ► t]q and Kl — > vr°,r/8 transitions, respectively): 

1 [3 (vo |^ s=1 | K L ) 3 GI + 2G 27 /3 _ 3 x x 

P l - 2^ 2(n°\H%r S = 1 \K L ) 2 Gs -G 27 2 x 6 - x x ' 1 j 

Obviously, the nonet symmetry limit, p = 1, amounts to discarding the dominant Q\ contribu- 
tion (i.e. Kl —* uu — > 77), and keeping only the Qq one. Therefore, in that limit, these models 
should under-estimate the Kl — * 77 rate, but this is not the case (on the contrary, they need 
p < 1 to reduce the Qq contribution). In other words, from the dominance of the Kl — ► uu — > 77 
transition, a correct phenomenological pole model should be proportional to p— 1 (like Eg. (1441) ). 
to a good approximation. 



3.5 Comments on K L — > 77* 

For the Kl — ► decays, a detailed slope analysis was performed in Ref. [25], including 

vector meson exchanges. Two remarks can be made from the insight gained here. 

First, quite generally, our two-step phenomenological procedure to deal with pole amplitudes 
in U (3) could have some implications also for VMD-type models. Indeed, working first with 
theoretical masses permits the identification of vanishing contributions, for which it is very 
dangerous to include only partially the higher order corrections (due only to the masses for 
example). As a rule, cancellations plague pole amplitudes and should be dealt with carefully. 
This is especially true in SU(3) ChPT, since it is unable to catch a transition through a uu quark 
pair at leading order. Only once the surviving terms are precisely identified, one is allowed to 
switch to physical masses without risks. 

Given the specific topologies studied in Ref. [25], it is the parameter fr™ one * that could appear 
most sensitive to the previous remark. It originates from the processes Kl — ► Tr°,rj,r]' —* 
7(p, u>, 4>) — > 77*, which vanish at lowest order in SU (3) ChPT, and proceed through G| + IG27 
in £7(3) ChPT. However, all the dependences on the weak transition cancel in the normalization 
with respect to Kl — > 77, and by onet as given in Ref. [25] is not modified. 

Second, it is to be noted that the global sign of our Kl — > 77 amplitude when Gg/Gg < 
is opposite to the one obtained using phenomenological pole models 4 . This is supported by 
the analysis of Ref. [25]. Yet, a precise constraint on G|/Gg from Kl — ► 77* seems at present 

4 These models in general imply a destructive interference between the r\ and 7/ poles, leaving the 7r° one as 
dominant (see the discussion in [18]). Then, one can see (independently of our conventions) that the sign of the 
7T° pole in Eq.lj2HJ, keeping only Gs, is the opposite of the 77 one of Ea. l|45|l when Gg/Gs < 0. 
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Figure 5: Pole and tadpole diagrams for the process Kg — » 7t°jj. 



difficult to obtain given the many unknown counterterms and VMD couplings. So, the purpose 
of the remaining sections is to search for independent constraints on G%/G 8 from other radiative 
K decays. 

4 The Ks —> 7r°77 decay 

The radiative process Ks — ► 7r°77 is the simplest mode one can think of to test the physical 
mechanism advocated for Kl — > 77, with in particular Gg/G 8 — —1/3. The corresponding 
branching fraction has been measured two years ago by the NA48 collaboration, though with 
large uncertainties [26]: 

B (K S - vr°7 7 ) e m X 7 P 7>220 MeV = (4.9 ± 1.8) x lO" 8 . (50) 

The decay proceeds through pole diagrams at leading order [27], with also a tadpole graph in 
the case of the weak mass operator [28] (Fig. 5). Unlike for Kl — ► 77, the amplitude starts 
at 0(p 4 ) in both J7(3) and SU(3) ChPT. It can be parametrized in terms of a pole function 
B(T 2 ): 

A v» [K S -» vr 77) = ^B(T 2 )e^k lp k 2a , (51) 
with T = k% + k 2 , the total momentum of the photon pair. This gives for the width: 

r(tf fl -7rVy) = r* X dzA^(l,z,rl)zZ\B(z)F , (52) 

y[2Tr) J 0.2 

with z = T 2 jM\, X(a,b,c) = a 2 + b 2 + c 2 -2ab-2ac-2bc, r n = M^/Mk and z max = (1-rv) 2 ~ 
0.53. The infrared cut-off is introduced to get rid of the Ks — ► background as in Eq.(|50|). 

4.1 Pole amplitude in U(3) ChPT 

A tree-level 17(3) computation gives: 

B(T 2 ) = (m 2 K - ml) [G 8 B 8 (T 2 ) + G S 8 B S 8 (T 2 ) + G 27 B 27 (T 2 )] (53) 
with no contribution from Q™ and 



(m 2 K -T 2 )(m 2 -ml + T 2 ) 2 5m? K - 2m\ - 3T 2 

(T 2 -m 2 )(T 2 -m^)(T 2 -ml) ' s[ ' 2{T 2 - m 2 )(T 2 - m 2 ^) 

2 ^2\(rp2 ^2-- 



2 = 3m' (T 2 - m 2 K ) + 2(m' K - mj)(T 2 - m 
27{ ' 3(T 2 - m 2 ){T 2 - m 2 ,) {T 2 - ml) 
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Figure 6: a) B (K$ — ► 7r°77), x 10 8 , as a function of m for Gg/Gg = —0.3,0, +0.3 (dashed, 
plain, dotted), b) Comparison between the 0(p 4 ) SU(3) result (dotted), idem plus the m^ 2 
corrections (dashed), and full 0(p 4 ) U(3) result (plain), with G| = for the three curves. 



Note that l m^^,r\,rf refer to the theoretical masses, defined in Eqs.(|36j) and (|37jl. It is only when 
these masses are used (i.e., when one works consistently at lowest order) that the weak mass 
term Q™ indeed drops out. Note also the identity 

|s 8 - + ^B 27 = , (55) 

originating from the fact that Q2 cannot contribute to any of the subprocesses K$ — > 7r°(7r°, rj, rf). 

The behavior of the rate as a function of mo is plotted in Fig. 6a. A significant enhancement 
is observed when mo decreases. This is essentially an effect of the r\ pole, that comes closer 
to the upper boundary of phase-space when mo diminishes (Fig. 3a). In fact, this pole could 
even enter the phase-space if mo was chosen as low as mo < 447 MeV, which corresponds to 
m v < rriK — m n from Eq. (|HH|) . Such a value is of course to be avoided. The estimation of 
the actual enhancement requires the phenomenological prescription proposed in Sec. 3. 2, and is 
deferred to Sec. 4. 3. 

On the other hand, it is clear that, Qg being allowed to contribute at leading order, it 
will dominate the decay via the pion pole (and, to a lesser extent, the 77 one). This explains 
the moderate sensitivity of B (Kg — > 7r°77) to G| exhibited in Fig. 6a. However, the decay 
Ks — > vr°77 could still be a useful probe of Gg/Gg, as will be discussed in Sec.4.3. 

4.2 Reduction to SU(3) ChPT 

How can the predicted enhancement for decreasing mo (see Fig. 6a) be understood from the 
point of view of the SU(3) chiral expansion? In order to answer that question, let us develop 
the three pole functions (|54j) in powers of 1/m^: 

r> (T2 , _ T*-ml 2(5ml - 2ml - 3T*)(T* - ,^-4, 

8[ >~ (T 2 - m2 8 ) (T 2 - ml) 3{T* - m^m 2 + [ ' ' 

d (T 2^ T 2 -m\ {hm? K - 2ml - 3T 2 ){2m\ + m\ - 3T 2 ) , n( - 4) 

27 [ >~ (r 2 - m 2 g ) (T2 - ml) 9(T 2 - m^m* + 1 ' ' 

(56) 
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with m 2 s fixed by the GMO relation. Though Eq . (|54j) is just the 0(p A ) £7(3) amplitude, the 
various terms in the above series correspond to increasing p 2 orders in SU(3). As can be 
checked numerically, these series are quite well-behaved (i.e., as expected from naive chiral 
counting). Still, significant effects can build up through phase-space integration, leading to the 
large enhancement observed in Fig. 6a. The behavior of the total rate as a function of mo when 
the mo series is truncated at a given order is displayed in Fig. 6b: 

- Retaining only the 0(m{j) terms, the 0(p 4 ) SU(3) result [27] is recovered (dotted 
line in Fig. 6b): 

B [Kg - A7)S2'° (P4) = 3 - 8 >< 10 " 8 ■ (57) 

- Terms of order m^ 2 introduce a class of 0(p 6 ) effects, leading in 1/N C , which corre- 
sponds exactly to the Lj, N%s, L 9 and N± counterterms once these are saturated by 
770 exchanges (see appendix, Eqs. (j96l97ll00|0 . They are seen to conspire collectively 
at the rate level (dashed line in Fig. 6b), already without G|, in order to reproduce 
the effect of the 77 pole as a function of mo discussed before. 

- Keeping the full £7(3) result comes to keeping some leading 1/N C contributions at 
all orders. As the plain line in Fig. 6b shows, the leading 1/N C terms of order 0(p n>G ) 
amount to only a small correction to the rate (at least as long as mo is not too small) . 

4.3 The Ks —>■ 7r°77 rate and its sensitivity to the penguin fraction 

One of the advantages of dealing with 770 effects in the £7(3) framework rather than through 
SU(3) local counterterms is that we do not have to face the problem of fixing the mo parameter. 
Indeed, as discussed in Sec. 3. 2, the requirement of freezing the pseudoscalar poles at their right 
places, thereby restoring the analytical properties of the amplitude, leads to a well-defined 
prescription. Since both Q\ and contribute at lowest order, Eq. (|3§l) can be used directly, 
and we get to 



3{2M 2 K — M 2 — T 1 ) 2M 2 K + M 2 — 3T 2 I c e c n sgcy 



SiMl-M 2 ) \T 2 -M 2 T 2 -M 2 J 
+ 1T7k \ ^P?. kl 7 ). ~ ^9 'TIT ) ' ( 58 ) 



Bi (T 2 ) 



8(M 2 K 


-M 2 )(T 2 - 


-M 2 ) 


1 ( 


SgC v 


c e c v * 




T 2 - M 2 


T 2 - M 2 


3 ( 


s e c v 


CgC v > 


V2 \ 


T 2 - M 2 


T 2 - M 2 



2 _ -3(2M 2 K -Ml-T 2 ) 2M 2 K + M 2 -3T 2 / c eCr] s 0Crj , 



8(Mf C — M 2 )(T 2 — M 2 ) 8(M 2 K -M 2 ) \T 2 - M 2 T 2 - M 2 J ' 

with c v = eg — 2\[2s$ and cy = sq + 2\^2cg, the mixing angle combinations for 77 and 77' — > 77, 
respectively (as a consistency check, note that Ea.(|55|l is preserved). 

The resulting differential rate for z > 0.2 is shown in Fig. 7a. Yet its shape is not much 
affected by neither 6p nor Gg/Gg, and extracting information on the latter is clearly beyond 
foreseeable experimental sensitivity. Looking at the low-energy end of the 77 spectrum would 
not be more helpful. Indeed, integrating over < z < 0.048 (m 77 < 108 MeV), the obtained 
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Figure 7: a) K s -> vr°77 normalized differential rate for 9 P = -20°, G s s /G s = -0.3,0, +0.3 
(dashed, plain, dotted), b) B (Ks — > 7r°77) , x 10 8 , as a function of Gg/Gs for #p = 
-15°, -20°, -25°. The star refers to Eq.jHJ). 



partial branching ratio varies from 0.48 x 10 8 to 0.44 x 10 8 for G|/Gs between —0.6 and 0.6, 
while the 577(3) prediction is 0.49 x 10~ 8 . 

Fortunately, for the total rate, the situation is better though, as said before, G| is not 
dominant. Fig. 7b summarizes the constraints one could get on G|/Gg from an experimental 
determination of the branching ratio. Note that merely fixing the sign of G$/Gg (inaccessible 
from Kl — > 77 alone) would already give valuable information on the low-energy realization of 
the effective weak Hamiltonian (^Q). 

Of course, Fig. 7b is subject to the theoretical uncertainties associated with any leading 
order ChPT computation. A precise extraction of Gg/Gs from Ks — > vr°77 would require an 
estimation of the NLO effects. While some unitary corrections are already included through 
the use of Gs as extracted from Ks — ► vr7r, there remains the separate class of 0(p 6 ) effects 
from vector mesons. These were seen to give sizeable contributions to Kl — » 7r°77 (see Ref. [29] 
for a recent review and list of references) and though the underlying dynamics is very different 
for Ks — ► 7r°77, working out their impact would be really called for. In particular, one should 
investigate if they affect either the differential rate (Fig. 7a) or the total rate for low 77 invariant 
mass, which are pretty much insensitive to G%/G% and 9p. In that case, these two observables 
would provide essential discriminating tools, able to single out the 0(p & ) corrections due to 
vector resonances. 

Finally, one can predict B (Ks — > 7r°77) from the range of Gg/Gs preferred by Ea.(|22|) (which 
englobes the one extracted from Kl — > 77, see Eq.(|47|)): 

Gl/G s = -0.38 ± 0.12 B (K s -► ^77) ^ = ( 48 ± °- 5 ) x 10 ~ 8 » ( 59 ) 

where the error, inferred from Fig. 7b, includes the theoretical uncertainties associated with the 
physical mass prescription (i.e., with a class of 0(p 6 ) effects). As expected, the effect of the 770 
meson increases the 0(p 4 ) 517(3) prediction (J57|). Note that both results are in agreement with 
the (still imprecise) experimental value (|50j) . 

5 The K + —> 7r + 77 decay 

Let us now turn to K + — > tt + jj, which has already received a lot of attention. Experimentally, 
the situation has been improved recently [30] and is expected to get even better in the near 
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future. On the theoretical side, the analysis is slightly more involved than for K$ — > 7r°77 as 
there are now both loop and pole contributions at 0(p 4 ) [31]: 



A i»> ( K + ^+ 77 ) = |L \ A (T 2 ) (T 2 ^ - 2^^) + B (T 2 ) e^h^a 



(60) 



with yl(T 2 ), the Tr^-fT 1 * 1 loop function, and B(T 2 ), the ir°,rj,rj' pole one. These two pieces do 
not interfere in the rate since they produce the two photons in different CP states: 



n 2 ]\/fb f Z max r 

r (JT+- vr+ 77 ) = ±£k / dz X^(l, z ,rl) z *U\A(z)\ 2 + \B(z 

16 (47T) JO. 2 L 



(61) 



where the cut in z is applied to remove the K + — > 7r + 7r° background. Though the loop contri- 
bution is larger than the pole one, this latter piece should not be neglected. Indeed, as we will 
see, it is very sensitive to the ratio Gg/Gg thanks to the suppression of the tt° pole compared 
to the 77,7/ ones (K + — > 7r + 7r° is a A/ = 3/2 transition). 

5.1 Loop amplitude in SU(3) ChPT 

For the A (z) amplitude, working with the enlarged symmetry U(3) would be superfluous as only 
charged particles circulate the loop. We will thus take up the standard 0(p 4 ) SU(3) result [31]: 



A(z) 



G 8 

°27 — ~7T~ 



1 r l 
1 + *- + 5l 7 

z z 

13r 2 J_ 13 
3z 3z 3 



) $ (z/rl) + (l 



627 = G 8 



G 27 (7rl 13 13 



3z 



+ 3z 3 



(62) 



with the three-point loop functions 



d>(a) 



1 arcsin 

a 

1 

1 + - 

a 



l-Jl-A a . 
log V , + Z7T 



4/a > 1 

< 4/a < 1 . 



The contributions of the counterterms are collected in the unknown constant c: 



128tH 



c = 



3 (L 9 + L10) + iVu - N 15 - 2N 18 + 



2C27 
3G 8 



(3 (L 9 + L w ) + A) 



(63) 



(64) 



where the Lj's refer to the basis of 0(p 4 ) strong operators of Gasser-Leutwyler [32] and the iVj's 
to the octet \AS\ = 1 one of Ecker-Kambor-Wyler [33]. The 27-plet counterterms, suppressed by 
the AI = 1/2 rule, are collectively denoted by Di [6,34]. The L«, N{ and D- L occur in separately 
finite combinations in c. Finally, note that the weak mass operator does not contribute. This 
has to be so since it could not have been absorbed in the weak counterterms iVi4 15 is- 
The above loop amplitude induces the following branching fraction for z > 0.2: 



B (K + vr+77) 



L,SU(3) t O(P 4 ) 
2>0.2 



= (5.77 + 1.64c + 0.29c 2 ) x 10~ 7 . 



(65) 



This expression is not much affected by the momentum cut, as can be inferred from the shape of 
the two-photon invariant mass spectrum which exhibits the characteristic peak above the it + it~ 
threshold. 
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Figure 8: Pole and tadpole diagrams for the process K + — > 7r + 77. 



At 0(p 6 ), unitarity corrections from K —* tttttt as well as vector meson effects must be 
included, and were analyzed in Ref. [35]. The former increase the rate by about 30 — 40% while 
the latter should be smaller. From the corrected spectrum and rate, but assuming negligible 
pole contributions, a fit to the E787 data [36], (here we quote only the branching fraction) 



with i-L 



M K \ J2 „W 2 



\ CX P 



^ - 7j 100MeV< ^ +<180MeV = (6.0 ± 1.5 ± 0.7) x HP 7 (66) 
2 A (l, r 2 , z) L/z , the tt + momentum, leads to the 0(1) positive value 



1.8 ±0.6. 



(67) 



5.2 Pole amplitude in U(3) ChPT 



The pole diagrams for K + — ► 7r + 77 are similar to those for K$ — > 7r°77, with in addition a 
five-particle contact interaction from the WZW term 1)25(1 in the case of the weak mass operator 
Qg 1 (Fig. 8). A straightforward leading order computation in U{3) ChPT gives: 



with 



B (T 



b* It 



B 27 (T 2 ) 



m 



I) [G 8 B 8 (T 2 ) + G* 8 B S 8 (T z 



G27B27 {T 2 )] , 



10m^ + 3m 2 w + 7T2 



2(T 2 -m2)(T 2 
wiQ(5m| — m\ 



Bi (T 



%,) 
- AT 2 ) + (15m 2 K 



-5m 2 K + 2m 2 + 3T 2 
" (T 2 -m2)(r 2 -m2) 

7m 2 -8T 2 )(m 2 -T 2 ) 



(68) 



(69) 



3(T 2 -m 2 )(T 2 -m 2 )(T 2 



Again, m^^K^r)' stand for the theoretical masses, which ensures the absence of effects when 
the two tadpole graphs are included. Unlike for K$ — ► vr°77, the three pole functions B$, B% and 
S27 are now independent as Q±, Q2 and Qq all contribute. Note that the pion pole disappears 
in B§, leaving a constant term, because of the T 2 — m 2 momentum dependence of the off-shell 
K + — * 7T + 7r° amplitude (purely AI = 3/2 once on-shell). 

The pole-induced branching fraction as a function of mo is depicted in Fig. 9a. It is not a 
monotonically decreasing function down to the SU(3) limit as in the Kg — ► 7r°77 case. On the 
contrary, it exhibits a strong suppression for mo between 1.0 and 1.5 GeV. Below 1 GeV, the ry 
pole begins to be felt significantly, hence the rate increases as in Fig. 6a. 
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Figure 9: a) B {K+ -> 7r+77) po ' es , x 10 7 , as a function of m for G%/G% = -0.3, 0, +0.3 (dashed, 
plain, dotted) and for G% alone (dash-dotted), b) Comparison between the 0(p ) SU(3) result 
(dotted), idem plus the TOq 2 corrections (dashed), and full C(p 4 ) U (3) result (plain), with G| = 
for the three curves. 



5.3 Reduction of the pole amplitude to SU(3) ChPT 

The SU(3) limit of the above result can be investigated expanding the pole functions in powers 
of 1/m 2 ,: 

B (T 2^ -1/2 4(5m 2 Y -2m 2 -3T 2 )(m 2 Y -r 2 ) + n(m -^ 
§( )_ 3(T2 - m^fml 

B ^ = + 5 t^K + ° K " 4) ' 

R f^N _ 5/4 1/12 2(5m| - 2m 2 - 3T 2 ){3m 2 K - m 2 - 2T 2 ) _ 4 

*27{1)- T 2_ m 2+ T 2_ m 2 8 9(T2-m 2 8 )2m 2 +Uim ° j ' 

(70) 

The behavior of the total rate as a function of mo when the rriQ series is truncated at a given 
order is displayed in Fig. 9b: 

- The 0(p ) SU(3) amplitude (dotted line in Fig. 9b) corresponds to the 0(m[]) 
terms [31], and gives 

B (K + - ^)Z U 2 3) '° ip4) = L17 >< 10_7 • ( 71 ) 

This result is to be compared to 0.51 x 10~' without the G27 piece. Such a large 
effect of the 27 operator (see Fig. 9a) has been overlooked previously, and arises from 
the pion pole, absent from the octet, that compensates for the AI =1/2 suppression 
over the whole phase-space. 

- The ttiq 2 corrections are the collective effect of the SU(3) counterterms L7, N13, 
and ivj 6 \ when saturated by the 770 (see appendix, Eas. 1)9619711001) ). The series seems 
not well-behaved in the case of Bg, again because the leading term is suppressed by 
the absence of the ir° pole factor. In particular, for tuq between 1 GeV and 1.5 GeV, 
the tjIq 2 and mg terms interfere destructively and even cancel each other out for a 
given T 2 inside the phase-space. This is the origin of the dip shown in Fig. 9a, as 
confirmed by the dashed line in Fig. 9b. Said differently, the AI = 1/2 enhanced B$ 
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amplitude, being essentially due to the 77s pole, is very sensitive to r\ — rf mixing 
effects. 

- It is a well-known fact that when the leading order is suppressed for some reason, one 
can expect sizeable effects from the NLO corrections. Since in the present case the 
NLO corrections are not suppressed by any mechanism (compare the ttt-q 2 correction 
in Eas.([70j) and (|55jll. the series is expected to behave correctly afterwards. The 
plain line in Fig. 9b shows that this is indeed the case: the full U(3) result is rather 
well reproduced by the tuq 2 corrections alone. 

Note, finally, that for uiq < 1 GeV, B% and -B27 begin to interfere destructively, instead of 
constructively in the SU(3) limit. 

5.4 Pole contribution to the rate 

Given the strong sensitivity of the pole amplitude to the Tj — rj' system, the prescription Eq, (|39j) 
has to be enforced before numerical evaluation. Expressions similar to those of Ea. (|58|) are then 

3 2M 2 K + Ml - 3T 2 / c c v , sjc^ \ 

4(M 2 K - Ml) + 4(M 2 - Ml) \T^- M* + T 2 - M 2 ) 

J_ ( S d C v _ C d C v , \ 

3_ / s c v c e c v > \ 

2^2 \T*-M* T^-M 2 ) ' 

5M| - 7 Ml + 2T 2 3M 2 - — M 2 — 2T 2 ( c e c n s eCrf \ 

4{M 2 K - Ml)(T 2 -Ml) + 4(M| - M%) \ T 2 - Af 2 + T 2 - M 2 J ' 

with c v = eg — 2\[2sq and c v > = sg + 2\/2cg, as in Sec. 4. 3 5 . The pole contribution to the rate 
as a function of G|/Gs is depicted in Fig.lOb for z > 0.2. 

In the positive sign alternative for GyGs, it could account for an increase of the total rate 
by more than 20%, and should be taken into account in the extraction of c, Ea. (|67|) . The shape 
of the differential rate could be affected too, but only mildly, as illustrated for a particular set 
of parameters in Fig. 10a. 

Conversely, for the preferred value G^/Gs — —1/3, the pole contribution is completely 
suppressed (about 10 -9 ). This is a pure numerical coincidence: the three contributions Qs,8s,27 
are of the same order, but can interfere destructively for G|/Gg < 0. No definite prediction can 
be made in this case since 0(p e ) corrections can no longer be neglected. Let us thus simply give 
an upper bound: 

G s 8 /G 8 = -0.38 ± 0.12 => B (K + -> 7r + 77)^ ( 2 ),0(p4) < 0.3 x 10~ 7 . (73) 

5 Note that, for both Ks — ► 7r°77 and K + — > 7r + 77, the two-mixing angle result is immediately obtained by 
subtituting (fc = sec (9o — 0s) F n /Ft) 

cety — > foCa(fss + 2V2f cs), sgc v / — » /gSo(/sSo + 2^2/o c 8 ) 
s e c v — > fos s (fsco - 2V2/0S8), cgc v — > facoif&ca - 2V2f s s ) 



obtained: 



B 8 (T ) 



bi (r ) 



B27 (T 2 ) 
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Figure 10: a) K + — > 7r + 77 normalized differential rate for 0(p 4 ) loop + pole (plain), loop alone 
(dashed) and pole alone (dotted) contributions, with c = ±2 and G^/Gg = +1/3, 9p = —20°. 
b) B(K+ -> 7r+77) poies , x 10 7 , as a function of Gf/G 8 for P = -15°, -20°, -25°. The star 
refers to Eq. ([23|) . 



This could go up to around 0.5 x 10 7 with C(p 6 ) effects but, in any case, no signal of pole 
contribution should show up experimentally, neither in the total nor in the differential rate. 

5.5 Analysis of the low energy end of the 77 spectrum 

Recently, the E949 Collaboration has obtained the following upper limit [30]: 

B (K+ - vr+77)p; +>213 MeV < 8 " 3 x 10 ~ 9 > ^ 

assuming a spectrum of the shape predicted by ChPT with unitarity corrections. This corre- 
sponds to a 77 invariant mass in the range < z < 0.0483, below the 7r° pole. The loop and 
pole leading order SU(3) predictions for this range are given by: 

B [K+ -> vr+77)t<o' 7 o ( 4 3 83 <!,(p4) = I - 062 ' + 0A5d + °- 09 ) x 10 ~ 9 ' ( 75 ) 
B (K+ -> vr^^^oS ^^ = L0 >< 10 " 9 ■ ( 76 ) 

For c of order 1, the loop contribution should be < 10 -9 . Including the unitarity corrections 
analyzed in Ref. [35], which are specially large for small z, it raises to 

B {K+ -> ^ + 77)S 3 83 ,im<a " 2e " = 6- 1 x 10_9 ' ( 77 ) 

A large error (presumably more than 50%) should be assigned to this number given the rather 
poor theoretical control over this small corner of phase-space. 

The pole contribution, including 7/0 effects, should also be taken into account. Restricting the 
phase-space integration to the present range, we find the situation shown in Fig. 11. Interestingly, 
it appears that given the rather large unitarity corrections, and the present experimental bound 
on the high 7r + momentum end of the spectrum, negative values of Gg/Gg are preferred. Together 
with the possibility of correlated study of the total and differential rates, this makes of K + — ► 
7r + 77 a promising mode to gain information on this ratio. 
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Figure 11: B(K+ -> 7r+77) po ' es for z < 0.0483, x 10 9 . Assuming non-negligible loop contribu- 
tions [35], the recent upper bound [30] hints towards negative values for Gg/Gg. The star refers 
to Eq. lj23|) . 



6 Other radiative modes 

In this section, we consider two other modes involving pole amplitudes. Our discussion will be 
very brief, because the sensitivity to the parameter of interest Gg/Gg will turn out to be either 
very suppressed like for Kl — ► 7r°7r°77, or buried among other unknown parameters like for 
Kl —> vr + 7r _ 7. 

6.1 The K L — > 7T 7r°77 decay 

This process is entirely due to pole diagrams at 0(p ) (see Fig. 12), and was analyzed in SU(3) 
ChPT in Ref. [37]. The U(3) amplitude is very cumbersome and will not be given explicitly. Its 
main features are (independently of the use of the physical mass prescription): 

- It receives contributions from Q$, Q| and Q27- As for Kl — ► 77, there is no tadpole 
graph but the contribution of still cancels out when the rj, rj theoretical masses 
are inserted. 

- It depends on both the T 2 (77 invariant mass) and Q 2 (tttt invariant mass) kine- 
matical variables, while the SU(3) amplitude depends only on T 2 . 

- As for Ks — > 7r°77, the G%,G\ and G27 pieces are not independent because the 
current x current Q2 operator does not contribute. 

At first sight, one could think that the possible double occurrence of the pseudoscalar poles 
would give a good sensitivity to Gr|. Unfortunately, this is not the case because the available 
phase-space for the 77 invariant mass is very much concentrated around the 7r° pole. Since 
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77, rj — > 7r 7r 7r° are both isospin violating, it is then again the 7r° pole that dominates the initial 
Kl — > -P transition. The effects of the 77, rf and Gg are thus completely suppressed. 

Numerically, we have checked that, indeed, dealing with the mode in U{3) ChPT with the 
physical mass prescription Eq.(|39j) changes the SU(3) result of Ref. [37] only at the percent level, 
no matter the 77 invariant mass cuts. In addition, differential rates with respect to either T 2 or 
Q 2 are only slightly modified. In conclusion, the effects of the 770 and G§ are at the percent level, 
beyond experimental sensitivity, and also beyond theoretical control since 0(p G ) corrections to 
the dominant tt° pole amplitudes are certainly bigger than a few percents. 

6.2 The K L — > 7r + 7r 7 decay 

This mode receives many types of contributions, so let us briefly describe its structure. The decay 
amplitude can be decomposed into electric and magnetic transitions. The electric piece starts 
at 0(p 2 ), and is completely determined at this order by the bremsstrahlung from Kl — > 7r + 7r~. 
It is therefore CP-violating and suppressed. At 0(p 4 ), there is also a loop contribution, which 
is quite small [38]. 

This has allowed experimental access to the magnetic transition [39,40]. This one starts at 
0(p A ) with contributions from local AS = 1 odd-parity operators (~ N29 + A31, see Ref. [38]) 
and pole diagrams (Fig. 13a). In addition, 0(p e ) effects were seen to be important from the 
experimental observation of a non-zero slope in the photon energy. A detailed study at that 
order is done in Ref. [41], including both vector meson exchanges and chiral loops, and we will 
rely on that work for conventions and numerics (as well as for a complete reference list). 

In the present work, we wish only to comment on the pole part of the magnetic amplitude 
(the F\ piece in Refs. [38,41], shown in Fig. 13a), which has been seen to be very sensitive to 
nonet symmetry breaking [19]. In our framework, this translates into a sensitivity to G|. 

As for Kl — > 77, the pole amplitude vanishes in SU{3) at 0(p i ) when enforcing the GMO 
relation. In U(3), it is simply obtained by substituting VJfj in Eq. (|3*2*|) by 

V+ _ = jJ- I I ie^pfplkl, (78) 

n n 7 2^/2vr 2 F 3 I J v p a 

and thus arises again only through the uu component of the weak vertex, i.e. through Q\. As 
for Kl — ► 77, NLO effects from Qq are assumed to behave according to usual chiral counting 
and thus are discarded. Note also that the analysis of Kl — ► 77 at 0(p e ) in SU(3) presented 
in the appendix would not be much modified for the present mode, and that its conclusions are 
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Figure 13: a) Anomalous pole contribution to the magnetic Kl — * tt^tt^j amplitude, b) 
Numerical correspondance between the parameter p defined in [41] and G^/Gs (see text). The 
star refers to Ea. ()2"3*)) . 



still valid. For these reasons, the phenomenological pole parametrization of F\ usually adopted 
is not appropriate and should be replaced by 

gg + |G 27 / 1 (c e -V2s e f ( Se + V2c e f \ 

G 8 \Ml -M2 + 3(Ml - M2) + 3(M| - Mj) ) ' U9j 

which has no contribution, is now dominated by the rj pole and is thus negative when 
G s 8 /G 8 < 0. 

Instead of repeating the numerical analysis of Ref. [41], we have chosen to simply give a 
dictionary relating their parameter p to the ratio Gt/Gg for fixed F\ (Fig. 13b). In other words, 
to each value of F\ corresponds a given p in Ref. [41] , and a given G| in Eq. 1)79)1 , which are then 
reported in Fig. 13b. It is important to understand that this is only a numerical equivalence, 
and does not correspond to the analytic equivalence of Eq.(jl9*)l (for example, G%/G% = —1/3 
implies p ~ 0.5 via Eg. 1)15)1 . but is accounted for by the smaller effective value p = —0.15 ± 0.07 
in Fig. 13b). As already explained in Sec. 3. 4, the usual phenomenological pole model is incorrect 
because not proportional to p — 1 (to a good approximation). 

Unfortunately, in Ref. [41], p was only marginally allowed into negative territory because 
this was believed to be in contradiction with Kl — ► 77. As we have seen, Kl — ► 77 points 
towards Gq/G$ — —1/3, and therefore negative values of p should be allowed. From the fits of 
Ref. [41], it seems that such values could accommodate the data, but additional work would be 
needed. In any case, it is clear that a precise extraction of Gg/Gs from Kl — * ir + ir~j would 
be quite intricate due to the presence of many counterterms and vector meson couplings whose 
estimations introduce some amount of model dependence. 

7 Compatibility with hadronic observables 

This final section concerns non-radiative observables. First, we will see how our understanding 
of Kl — > 77 can help estimating the 0(G 2 F ) Kl — K$ mass difference generated by pseudoscalar 
pole contributions. Then, we will turn to the direct CP-violating parameter e'/s, and analyze 
how the information gained on the AI = 1/2 rule fits in the usual theoretical analysis. 
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7.1 Pole contributions to AMls 

Experimentally, the Kl — K$ mass difference, AMls — M (Kl) — M (Ks), is quite well-known 
[24] 

AMf s p = (3.483 ± 0.006) x 10~ 12 MeV . (80) 

Unfortunately, theoretical control on AMls ~ 2 Re (K° \Hw\ K ) is a long-standing issue since 
both short-distance and long-distance contributions are present [42] 

AMls = AMff + AMff . (81) 

The first piece, from VF-box diagrams with virtual c, t quarks, accounts for the bulk of the 
experimental value, AM^/AM e L J = (86 ± 26) % [43] (see also [44]). The second one arises 
from VF-box diagrams with low-virtuality u quarks, hence has to be dealt with light meson 
exchanges. In the present work, we will concentrate on the long-distance, non-local contributions 
from double insertion of \AS\ = 1 transitions (see for example [42,45]). 

In SU(3) ChPT, only pole and tadpole diagrams (Fig. 14a) occur at 0(p 2 ), giving rise to the 
0(G 2 F ) mass difference: 

AMlf SU ± S) 2F*m* K (G 27 -G 8 + G™) 2 ( 2 1 2 + 1/3 2 ) . (82) 

Though both M (Kl) and M (Ks) are renormalized by a term in (G™) 2 from the tadpole 
diagrams (Fig. 14a) 6 , the mass difference is independent of Gg 1 as can be seen by enforcing the 
GMO relation (|29|) . Still, everything cancels along with it, and AM^I? exactly vanishes. 

Pseudoscalar poles thus begin to contribute at 0(p A ) in SU(3), through loop corrections to 
the masses and weak vertices, along with genuine two-pseudoscalar loops, K — > PP — > K, from 
which they cannot be disentangled. Overall, these 0(p 4: ) loops require some unknown 0(p A ) 
local counterterms, including \AS\ = 2 ones, to cancel their divergence (see the comment at the 
end of the appendix). Alternatively as they correspond to the low-energy tail of the H^-box 
diagrams, an approximate matching with short-distance can be implemented [46]. 

Performing the same analysis in U(3) has the immediate advantage that the transition sd — > 
uu — ► sd is caught at leading order (Fig. 14b). Adapting Ea. H32|) . with Vi^ replaced by V we . a ki 
we reach 

AMls ~ m 2 - 3m\ + 3ml {° 8 + 3° 27 ) ~ m 2 - 3m\ + 3m\ { ° WXl) ' (83) 

Though Qs and Q™ cancel again upon enforcing the theoretical r/o — rjs mass matrix (|31[l. the 
non-zero contribution from (Qi) 2 survives (neglecting Q3 and Q5). As for Kl — ► 77, Q2 and 
<56 cannot contribute at leading order (such a disappearance of Qq has already been noticed in 
Ref. [20]). The physical mass prescription is then used to restore correct analytical properties, 
but only for the (Q1) 2 contribution 



AMlf = ^M\ ( G% + 2 -G 27 ) 2 ( a 1 , + {C ; ~ ^ Sd) l + { * e ± ^ C9 ^ I . (84) 
J ' s K 1 8 3 7 I M 2 - - M^ 3(M| r - M 2 ) 3(Af£ 1 fi 1 1 




6 In the usual language (see Ref. [5]), this is a manifestation of the necessary realignment of the vacuum brought 
by the tadpole operator Qg 1 . As a result, at 0(G%), the physical Kl.s masses are free parameters, and only the 
mass difference is calculable. 
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Figure 14: a) Pole and tadpole diagrams for AMff. b) Long-dist ance uu contribution, c) 
Fraction of pole contribution to AM^c as a function of Gg/G$ for Op = —15°, —20°, —25°. The 
star refers to Eq. (|23]) . 



Once again, the rj contribution dominates such that poles give a negative contribution to AM^s, 
which grows (quasi) quadratically with Gg/Gs (Fig. 14c). Note also the similarity of Eq.(|84|) with 
Eq.(J7nj), stemming from the proportionality of the vertices V^ +7r _^ ((75)) and V wea k when 
x 6 -> 0. 

To reach Eq.(JHIJ), we have discarded x§ because it does not occur at lowest order, i.e. in 
Eq. (|83]) . As explained in Sec. 3, in such cases, the physical mass prescription would generate 
some contributions corresponding to the inclusion of only a small class of higher order effects, 
and pole amplitudes being plagued by cancellations, this is a very dangerous procedure. In the 
case of the Kl — > 77 (and Kl — > tt + tt~^) radiative mode, we went on by arguing that Qq is in 
fact suppressed at higher order because the quantum numbers of the initial Kl together with the 
electromagnetic WZW vertex project out non-nn transitions, and kept only the Q\ contribution 
in our comparison with experiment. 

This last step cannot be extended to AMls since there is no such projection from CP- 
symmetry and charge conservation for the K° — K° transition. In particular, the (Qe) 2 contri- 
bution at C(p 4 ) is certainly not small now since it contains part of the K — > PP — » K loops, 
i.e. the low-energy tail of the box diagram depicted in Ref. [46]. These loops give a positive 
contribution to AM is, and though their estimates vary greatly, they are generally of a few tens 
of percent of AM^ P . 

Eq. (|84j) is therefore the correct leading order estimation for the long-distance (Q1) 2 piece 
only. To compare it with experiment, let us take, with respect to AM^ P , the conservative upper 
bound of +50% for the K — > PP — > K loops, assuming they saturate the (Qe) 2 contribution, and 
+150% for the short-distance piece. We can then safely infer from Fig. 14c that pole contributions 
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should be no less than —100%, i.e. 



\G S 8 /G 8 \ < 0.6 . (85) 

Furthermore, for our preferred value Gg/Gg ~ —1/3, poles contribute for about —30% such that 
the long-distance contributions partially cancel each other. 

7.2 Strong penguin contribution to e'/e 

The direct CP-violation observable e'/e is related to the imaginary part of the K — > irir isospin 
amplitudes, and additional theoretical inputs are necessary. We do not intend to give a full 
account of e'/e theory (for recent theoretical updates and references, see [47]), but would like 
to see if the large penguin contribution to the AI = 1/2 rule found at the hadronic scale is 
compatible with the small e'/e observed [48-51] 



(1.67 ± 0.16) x 10^. (86) 



cxp 



The general formula for this quantity is, in terms of the isospin amplitudes defined in Ea. (|16|) . 



e 



V2\ 



Im A 2 Im Aq 



Re A 2 Re A 



with $ = -8 + ^ « . (87) 
4 



We will concentrate on the strong QCD penguin contribution, and thus discard isospin breaking 
effects due to 7r° — 77W and electroweak penguins. Keeping then only the ImAo/KeAo part, we 
get in terms of matrix elements of the dominant Qi, Q 2 and Qq four-quark operators 

e'\ u lm\ t 3y 6 (m) {Qq (a)} 



y/2 \e\ Re\ u -z x (fJ,) (Qi (n)) + 2z 2 (fi) {Q 2 (m))o + 3z 6 (m) (Qe (m))o ' 

with (Qi) Q = (irir (I = 0) \Qi\ K) and \ = V ld V*. 

At this level, the usual treatment consists in taking the experimental value for the denom- 
inator (i.e. Re ^4o)) an d trying to modelize the Qq (fi) matrix element in the numerator at a 
scale for which i/q (/i) is calculable. In this respect, the penguin fraction Tp does not help much 
since running it up to /i > 1 GeV throughout the non-perturbative regime is beyond our reach. 
Alternatively, to make use of the extra theoretical and phenomenological information Tp — 2/3, 
see Fig. 4, let us rewrite Ea. (|88|) in a way independent of the hadronic matrix elements [52] 

_ u Im At ^ y 6 (n hadr ) 



ej \/2\e\ ReA u z & (^ had r) 

with Hhadr the typical scale of ChPT. The difficulty now is shifted from getting up to (Qq (//)) 
in the perturbative regime, to getting down to the ratio ye (/x) / zq (/x) in the non-perturbative 
regime. 

From the experimental constraints on the CKM factors Aj, it follows 



^ ) z 6 {fJ'hadr . 



Allowing for at most a 50% reduction due to isospin breaking effects, we obtain then 

\<Tp V -^^<\. (91) 
5 z e {fi had r) 5 
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Figure 15: The running of the ratio y§{ij) / zq{[i) at leading and next-to-leading logarithmic 
approximation, in the MS / NDR scheme, for as(Mz) = 0.118, m c = 1.3 GeV (based on [2]). 



For Tp ~ 2/3, this range is by no means unrealistic. As shown in Fig.15, the ratio y% (//) / zq (//) 
exhibits a scale dependence (this is obvious since zq (m c ) = at LO), which quickly decreases 
with [i [2] . This behavior arises because at both LO and NLO in the NDR scheme, y§ (/i) — zq (fj,) 
is scale independent to better than 10% for /x even as low as 0.6 GeV where clearly perturbation 
theory should no longer be trusted. It is therefore possible that though y§ (fJ-hadr) and zq (nhadr) 
are completely beyond our reach, their ratio could be controlled. If, for once in the long story of 
e'/e, nature was kind enough to somehow protect the specific ratio ye (fi) /zq (fi), an alternative 
theoretical strategy would be revived. This remains to be seen though all we can say (on the 
basis of Fig.15) is that a value smaller than one for this ratio is certainly not ruled out, and 
would require the bound 

G s 8 /G 8 > -0.8, (92) 

which is not much more constraining than what would have been obtained from the QCD- 
inspired analysis of Sec. 2. 3 (see Fig.l), i.e. Gg/Gg > —1. 

8 Conclusion 

A number of radiative K decays involving pseudoscalar pole diagrams have been investigated 
in the context of large- N c ChPT. Emphasis has been laid on the AS" = 1 weak operator Q|, 
peculiar to the U(3) framework, that holds the key to a better understanding of the underlying 
flavor-changing mechanisms and, in particular, to a phenomenological extraction of the penguin 
contribution to the AI = 1/2 rule. Let us now summarize our results, divided into three 
categories: 

Phenomenological constraints on the weak coupling G| from radiative K decays 

- The Kl — > 77 decay turns out to be mainly driven by the nonet-symmetry breaking 
operators Q| and Q27 (i-e., Qi), and dominated by the 77 pole. The contribution of Qs 
(i.e., Qq) is suppressed by large cancellations. Experimental data then imply G^/Gs — 
±1/3. The analysis of Kl — > 77* of Ref. [25] favors the negative solution. This sign is 
important for the interference between the short-distance and dispersive 77 contributions 
to K L fj,+fj,- [18,53]. 
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- The Ks — > 7r°77 rate is enhanced by r/o effects for most values of G%/Gg. In addition, the 
fact that neither the shape of the m 77 spectrum nor the rate for low m 77 are affected by 
G§ could be exploited to disentangle the NLO effects (like vector resonances) from the ijq 
ones, and to achieve a clean extraction of G\/Gg. 

- In the case of K + — > 7r + 77, the pole contribution to the total rate is negligible with 
respect to loops for Gg/Gg < 0, but can be as large as 20% for Gg/Gg > (poles are then 
needed to extract c). This situation arises from the fact that the Gg, Gg and G27 pole 
contributions are of the same order. At the low energy end of the 77 spectrum, the recent 
experimental upper bound [30] points towards Gg/Gg < (still, theoretical uncertainties 
are large in this small corner of phase-space) . 

- The effects of 77, rf and Gg on Kl — > 7r°7r°77 have proved negligible. 

- Relying on Ref. [41], the Kl — ► 7r + 7r~7 data are found compatible with a large range 
of Gg/Gg values. A precise extraction seems beyond reach given the many unknown 
phenomenological parameters. 

- The pole contributions to AM^g were also briefly discussed. As for Kl — > 77, they 
arise from Q| and Q27 (i.e., Q\). Yet, this time, non-negligible Qg (i.e., Qq) effects are 
expected at NLO from the presence of K — > PP — > A" loops. The poles result in a negative 
contribution to AMls, while loop and short-distance W-box. diagram contributions are 
positive. Conservative bounds on the latter two lead to |G|/Gs| < 0.6. For Gg/Gg ~ —1/3, 
large cancellations are expected between the loop and pole long-distance contributions. 

QCD-inspired operator basis and the penguin contribution to the AI = 1/2 rule 

- The three weak operators (Qg,Qg,Q27) of 17(3) ChPT can be related to the low-energy 
realizations of the current-current and penguin operators (Qi, Q2, Qg)- QCD-inspired the- 
oretical expectations then lead to the range (Gg/Gg) th = —0.38 ±0.12, compatible with all 
the phenomenological constraints given above. The non-perturbative evolution of current- 
current operators is therefore rather smooth, while the Qq penguin is significantly en- 
hanced, and responsible for about 2/3 of the AI = 1/2 rule at the hadronic scale. This 
large penguin fraction is compatible with the small e'/e observed, though a significant 
cancellation with isospin breaking effects might be welcome. 

- For Kl — > 77, Kl — > it + it~j and AM^I?, the well-known vanishing of pole contributions 
at lowest order in SU{3) ChPT and the subsequent pathological sensitivity to NLO effects 
are explained. Indeed, the transitions Kl — ► uu — > ^/tt + -k~^/ Kl, i.e. through Q±, 
cannot be caught in SU (3) at LO simply because there are not enough independent weak 
operators. This does not occur in U (3), where the uu leading order effect can be identified. 
Contributions from NLO (in particular from the dd, ss transitions driven by Qq) should 
then behave according to chiral counting, except for AM^g where K — > PP — > K loops 
are present. 
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Theoretical progress in the treatment of pole amplitudes 

- We have proposed a two-step procedure to deal with pole amplitudes: first, the iden- 
tification of vanishing contributions through the enforcement of theoretical masses (i.e., 
working consistently at a given order) , then the restoration of correct analytical properties 
by setting the poles at their right places for the remaining contributions only. 

- The weak mass term has been shown to disappear in all the modes considered, in both 
SU (3) and U (3) ChPT, as long as one is working consistently at a given order (see above). 
This is non-trivial since pole amplitudes involve both off-shell weak transitions [28] and 
the WZW action. 

- The connection between U(3) and SU(3) ChPT has been analyzed in details for all modes 
and, though sometimes large, 7]o effects are not incompatible with naive SU(3) ChPT 
power counting. Indeed, they can be reproduced, to a good approximation, with four 
SU (3) NLO counterterms saturated by 770 exchanges. For Kl — ► 77, a complete analysis 
at 0(p 6 ) has been presented, showing in particular that decay constant corrections do not 
occur, and that saturating the Kl — ► 77 amplitude with Q\ amounts essentially to assume 
that a combination of scalar-dominated counterterms is small. 

In conclusion, a consistent picture seems to emerge from our analysis. In particular, naive 
expectations from QCD are already well-supported by phenomenological constraints. Additional 
experimental information from radiative K decays are eagerly awaited, as they could give further 
insight into the interplay between strong and weak interactions at low energies. 
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Appendix: K L -> 77 in SU{3) ChPT at 0{p 6 ) 

In the text, we argued that corrections to the weak vertices Kl — ► ^ ,^8i to the 7r°,r/8 prop- 
agators and to the vr ,r/g — > 77 vertices (in particular through F — > F^,F V ^) should either 
reconstruct the dominant sd(ds) — ► uu — > 77 transition, or vanish (to a large extent). This is a 
strong statement since it implies a very specific interplay between the weak and strong sectors of 
SU(3) ChPT, and between 0(p 4 ) and 0(p 6 ) counterterms. It is the purpose of this appendix to 
check this assertion by performing the full calculation of Kl — > 77 in SU{3) ChPT at C(p 6 ). In 
particular, we will detail all the cancellations occurring between the various corrections depicted 
in Fig. 16. 
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Full 0(p 6 ) amplitude 

From the diagrams of Fig. 16, summing over the tt° and r]$ poles, we obtain: 



(93) 



with 
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(94) 



Aj = m 2 (D £ - logm 2 //x 2 + 1) and A 



7 + log47r. We have used the standard basis of 



Gasser-Leutwyler [32] for the 0(p 4 ) strong counterterms (A) and the Ecker-Kambor-Wyler [33] 
one for the 0(p 4 ) octet AS = 1 counterterms (iVj). In the third diagram, the CT (countert- 
erm) contribution originates from the odd-parity 0(p e ) strong Lagrangian, which is written in 
the Ebertshauser-Fearing-Scherer [54] or Bijnens-Girlanda-Talavera [55] basis (up to an overall 
normalization): 
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(96) 



with x± = U\ i xU^ , f± = F£ v ± UF^U^ (in our case x = r M with M the light quark 
mass matrix, Fj* u = F^ = —eQF^ v with F^ v the QED field strength tensor and Q the light 
quark charge matrix). In the last diagram, the CT contribution originates from the odd-parity 
0(p e ) weak Lagrangian. To our knowledge, no operator basis has been set up for it yet. For 
the purpose of the present analysis, only one of the direct Kl — * 77 couplings will be needed 
explicitly: 
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2^1 



All the others are collectively denoted by iV^ (octet part). Note that the C(p 4 ) CTs Lj and Ni 
alone suffice to absorb the Gs part of the loop divergence, Lf ] and Af } being thus separately 
finite. Finally, also the 27-plet CTs, though needed for renormalization, will not be specified 
explicitly, and are collectively denoted by A [6, 34] and D\ . 

Note also that, to reach Ea. (|94j) . the initial Kl mass Mj^ has been expressed back in terms of 
the bare mass m 2 K (see Ref. [32]). Alternatively, the weak rotation can be performed to discard 
the K L -> 7r°,77 8 O (p 2 ) vertices [27,31]. 
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Figure 16: Kl —>■ 77 at O (p e ). Grey vertices stand for local counterterms and meson loops. 



Discussion of the counterterms 



A first encouraging observation is that many counterterms drop out when summing over the 
7T° and 778 poles (^5^,10,11 and £4,5,6) or taking the external particles on-shell (like ^29,-^31), 
while the final combination of chiral logs is very simple and very small (it vanishes for \x ~ 550 
MeV). In particular, the a priori significant correction due to SU(3) breaking effects in the decay 
constants (i.e., £4,5) drops out, the impact of F n s/F n ^ 1 being thus at least of 0(p 8 ). But this 
is not the end of the story: there are still many interplays between the counterterms, though it 
is necessary to go beyond the strict SU(3) framework to get hold of them. 

As is well-known [32], L7 is well-described by a tree-level 770 exchange, which gives it an 
abnormal counting L7 ~ N 2 (still, see the discussion in [56]). In the same way, ./V13, 
and (which also have zero anomalous dimensions) can be saturated by the Vo resonance, 
surviving then as well in the large- N c limit. The explicit reductions of the U(3) weak operators 
to 5*7(3) read (fj = ^/2/3rj /F) 



Qs = 4 [<A 6 L^> - F 2 (XqL )j ) (D"fjo)] , 
* - -6F 2 (X 6 L lM ) (D"fja) , 

F 4 [cos (770) (A 6 x+) + i sin (770) (A 6 X-)] 
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Q 



(98) 



where the left-hand side £ 17(3) and the right-hand side G SU(3) (Q27 does not couple to the 
singlet). Similarly, the WZW Z7 (3) Lagrangian is reduced as [57] 
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A straightforward computation then gives 
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(100) 



after partial integration and use of the SU(3) equations of motion. Note that other counterterms 
(like N- 



21, 



(6) L (6) 



11 



could be saturated by the 770, but none contributing for on-shell photons. 
Besides, we will set /4 6) = 0, as suggested by the success of the LO description of it —* 77. 
There is then only the counterterm for r\ — ► 77 at 0(p 6 ) (chiral loops and wavefunction 
renormalization amount to F — > F„s, see Ref. [58]). Using its 770-saturated value with mo ~ 
850 MeV, we obtain r (77 — ► 77) ~ 0.4 keV, instead of 0.1 keV without Experimentally, 
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r (rj — > 77) = (0.510 ± 0.026) keV [24], which means that indeed accounts for the bulk of 
NLO effects. This is simply a manifestation of the fact that U(3) at LO reproduces rj, rf — * 77 
reasonably well without extra operators like Eq. 1)26)1 . 

Inserting the values (J100|) in Eq. , we find that the Gg and G™ pieces of the r/o-saturated 

l + §< 



CTs indeed cancel out, leaving only a term proportional to G| + IG27 
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3 / uin 
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(101) 



Note that the use of the t7(3) lowest order constraint F = F n = F v q does not affect the observed 
cancellation since any deviation of F^q/F^ from unity can be accounted for by varying TOq 
(saturation from tree-level exchanges). 

The first term is of course exactly the U(3) result in the limit tuq — > 00 and corresponds to the 
Q\ contribution, i.e. to the transition sd(ds) — > uu — » 77. The G27 piece (Di,D^) is suppressed 
by the AI = 1/2 rule and can be discarded, along with the small chiral logs. The remaining 
0(p 4 ) counterterms Nq,N\2 and Lg can in principle give a large contribution to Kl — > 77 
because of their large numerical coefficients. Still, they all arise from scalar exchanges [33, 59] 
and it is reasonable to assume that they will behave similarly to the counterterms saturated by 
pseudoscalar exchanges, i.e. cancel among themselves to a large extent, or act as a correction to 
the leading Q\ piece. As a further clue, note that there must be some interplay between L$ and 
A^6,i2 in order to absorb the contribution of Gg" in the iVj's. Finally, the fate of the remaining 

0(p 6 ) local Kl — > 77 counterterms (N^) is not clear, but it is at least possible that they are 
small. 

In conclusion, we have pin-pointed the combination of counterterms on which we make an 
assumption when we saturate the Kl — > 77 rate by the Q\ contribution. We hope that some 
further work will be able to get at least an upper bound on this combination, so that the precise 
measurements of Kl — * 77 (|27[1 will reliably fix G\/G%. For now, as explained in the text, there 
are other theoretical and phenomenological indications that Gg/Gg has a value around —1/3, 
which then unambiguously implies that the remaining CTs in Eq. (|101|) must at least partially 
cancel among themselves. 



Comments on AMls 

Finally, let us make a few comments on the similar computation for AM^g. As explained in 
the text, in SU(3) at 0(p 4 ), there are both corrections to the pole amplitudes, similar to the 
Kl — > 77 ones given in Eq. 1)94)1 . and pseudoscalar loops like K — * PP — > K. Large cancellations 
among corrections are again observed, in particular among those for the decay constants and 
masses. It should be noted also that though the terms in G^Gg^ and (G™) 2 are non-local for 
M(Kl) and M{K$) (i.e., contain some chiral logs), they appear only locally for AM^g (i.e., 
they are proportional to 0(p 4 ) counterterms and can thus be absorbed in the iVj's). As at C(p 2 ), 
only the mass difference is unambiguous in Chiral Perturbation Theory (see Sec. 7.1). 

Concerning the reconstruction of the Q\ contribution, the role of and is taken up 
by the AS = 2 counterterm (XqX-) 2 , saturated by the r/o- Again, L7, ./V13 and this AS = 2 
counterterm appear initially with a complicated coefficient involving G§, but they conspire to 
reconstruct the sd — > uu — > sd transition, i.e. a term with the coefficient (Gg + 2G27/3) 2 . At this 
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stage, an expression similar to Eq. lJlOljl is reached, but with in addition some K — > PP — > K 
loop functions (and divergences), and a correspondingly different combination of the remaining 
counterterms. 

In conclusion, because of the presence of K —* PP — > A' loops and AS = 2 counterterms, it 
does not appear possible to get extra information from AM^J? on the counterterms remaining 
for Kl — > 77 in Eq. ljlOljL Still, the interplays among corrections leading to the reconstruction 
of the Q\ contribution are found to occur also for AM is- 
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